Microwave Photon Detector in Circuit QED 
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In this work we design a metamaterial composed of discrete superconducting elements that implements a 
high-efficiency microwave photon detector. Our design consists of a microwave guide coupled to an array of 
metastable quantum circuits, whose internal states are irreversibly changed due to the absorption of photons. 
This proposal can be widely applied to different physical systems and can be generalized to implement a mi- 
crowave photon counter. 



Quantum optical photodetection fll] has occupied a cen- 
tral role in understanding radiation-matter interactions. It 
has also contributed to the development of atomic physics 
and quantum optics, with applications to metrology, spec- 
troscopy, and quantum information processing The quan- 
tum microwave regime, originally explored using cavities and 
atoms ill 01, is seeing a novel boost with the generation of 
nonclassical propagatine fields jstl in circuit quantum elec- 
trodynamics (QED) (0. 0. Htl- In the last years we have wit- 
nessed a tremendous development of the field of quantum 
circuits [13, Eitl- These devices are built, among other 
things, from superconducting elements, Josephson junctions. 
Cooper-pair boxes lll2ll . SQUID's, microwave guides and cav- 
ities iHSSl, all of them cooled down to the quantum degen- 
erate regime. Among numerous applications we may high- 
light the creation of artificial atoms or circuits with discrete 
quantum energy levels, and quantized charge 11211 . flux lll3ll 
or phase ll4ll degrees of freedom. These circuits find applica- 
tions not only as quantum bits for quantum information pro- 
cessing as charge flux iflil [19I1 or phase qubits 




I 20[ 21, 221, but also in the linear and nolinear manipulation 



of quantum microwave fields. In particular we remark the ex- 
change of s ingle p hotons between superconducting qubits and 
resonators tzl |23l |24|, |25|l, the first theoretical efforts for de- 



tecting incoming photons 1260, the generation of propagating 
single photons [5] and the nonlinear effects that arise from the 
presence of a qubit in a resonator l27ll28ll . 



While the previous developments represent a successful 
marriage between quantum optics and mesoscopic physics, 
this promising field suffers from the absence of photodetec- 
tors. The existence of such devices in the optical regime al- 
lows a sophisticated analysis and manipulation of the radia- 
tion field which is crucial for quantum information process- 
ing and communication. This includes Bell inequality ex- 
periments, all optical and measurement-based quantum com- 
puting, quantum homodyne tomography, and most important 
quantum communication and cryptography 101 ■ 

There are several challenges to design a microwave pho- 
todetector in circuit QED: i) Available cryogenic linear am- 
plifiers are unable to resolve the few photon regime, ii) 
Free-space cross-section between microwave fields and mat- 
ter qubits are known to be small, iii) The use of cavities to 



FIG. 1: Detector design. Sketch showing a one-dimensional waveg- 
uide (gray) coupled to a set of three-level absorbers (blue squares) in 
arbitrary positions. The microwave field (dashed line) excites coher- 
ently the state |0) to the upper state which decays onto a long- 
lived stable state. An analogous setup uses current-biased Joseph- 
son Junctions (CBJJ), in which a washboard potential confines two 
metastable states that can decay into a continuum of current states. 



enhance the coupling introduces additional problems, such as 
the frequency mode matching and the compromise between 
high-Q and high reflectivity, iv) The impossibility of perform- 
ing continuous measurement without backaction 12611 . which 
leads to the problem of synchronizing the detection process 
with the arrival of the measured field. We observe that there 
are related theoretical and experimental proposals that may 
solve some but not all of these problems. Most of them are 
oriented towards unitary or coherent evolution of the photon- 
qubit system, aided by nonlinear dispersive 124 12711 or bifur- 
cation effects ll29ll . But so far, to our knowledge, there has 
been no implementation of a microwave photon counter or 
even a simple microwave photon detector 

Our goal in this work is to design a metamaterial that per- 
forms single-shot microwave photodetection via iiTeversible 
absorption of photons llsoll . It resembles a photographic film: 
when a photon enters the device, there is a large chance that 
it will be captured, leaving the system in a stable and meso- 
scopically distinguishable state that can be observed a pos- 
teriori. More precisely, we propose a general setup based 
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FIG. 2: (Color online) Detection efficiency when absorbers are on 
resonance (real 7). (Top) Absorption probability vs. effective decay 
rate 7 = Vvg/V'^ in dimensionless units for a setup with A'^ = 1, 2, 4 
and 8 qubits (black, green, blue, red) either in cluster (dashed) or ar- 
ray (solid) regime. The error bars account for random deviations in 
the individual absorber properties, 7; , of up to 40%. When absorbers 
are close together, the efficiency is limited to 50%, while in the other 
case there is no upper limit. (Bottom) Detection efficiency vs. sepa- 
ration d in a periodically distributed array of absorbers. 



on a one-dimensional waveguide that passes by a set of pho- 
ton absorbers [Fig. [T|. We neither rely on any kind of cav- 
ity device nor aim at reaching a matter-field strong-coupling 
regime. The absorbers along the waveguide may be built with 
bistable quantum circuits similar to the ones used for imple- 
menting qubits. These circuits should be able to capture a 
photon and decay from an initial state |0) into a stable state 
\g) . In analogy to the photographic film, these irreversible 
events constitute the measurement process itself. The final 
step consists on counting the number of activated absorbers 
which is related to the number of photons in the Fock state 
which was detected by the quantum measurement. The count- 
ing should be done after the absorption process has finished, 
thereby avoiding any disturbing backaction on the incoming 
signal. Remark that the detection process is in this sense pas- 



sive and does not require any kind of control. 

It is possible to develop a simple model for the detection 
based on the previous rather general requirements. We model 
the absorbing elements as three-level systems, with an internal 
frequency uj and a decay rate F from the excited state |1) to 
the final stable state \g), see Fig. [T] Mathematically, this is 
described with a master equation for the density matrix p of 
the absorbers and the waveguide 



dt 



P - 



[H, p] + Cp. 



(1) 



The Hamiltonian contains terms for the absorbers and the ra- 
diation fields, ipi and tp^, propagating left and right with group 
velocity Vg, see Methods. The interaction between both is 
modeled using a delta-potential of strength V 

H = ^huj\i)^(i\+ihvg dx v/a^Vi - V'J'S^jVr 

i 

+ J2 f dxVS{x - X^) [(^1 + ^r) (0| + H.C.] ,(2) 

where Xi and |0) , |1) denote the position and the states of the 
i-th absorber The Liouvillian £ = has the standard 

decay terms for each of the absorbing qubits 



and it is proportional to the decay rate F. The solutions of 
this master equation can be found using an equivalent non- 
Hermitian Hamiltonian = H-i Yjj r/2 (1| , that rules 
the dynamics of the populations in |0) and |1) . The norm of 
the wavefunction is not preserved by this equation, but pre- 
cisely the decrease of the norm is the probability that one or 
more elements have absorbed a photon. 

The simplest scenario that we can consider is a single 
absorber coupled to the microwave guide, a problem that 
has analytical solutions for any pulse shape. In the limit 
of long wavepackets it becomes more convenient to analyze 
the scattering states of H . For a single absorbing element, 
these states are characterized by the intensity of the incom- 
ing beam, which we take as unity, and the intensity of the 
reflected and trasmitted beams, |rp and The associated 
complex amplitudes are related by the scattering matrix T, as 
in (t, 0)* = r(l, r)*. The single-absorber transfer matrix 



T 



1-1/7 -1/7 
1/7 1 + 1/7 



(4) 



is a function of a single complex dimensionless parameter 

7 = (r - i5)vg/V\ 



(5) 



which relates the properties of the circuit: the group velocity 
in the waveguide, Ug, the strength of the coupling between 
the absorbers and the waveguide, V, and the detuning of the 
photons from a characteristic frequency of the absorbers, 5 = 



3 



Lj ~ Lj^. The single-photon detection efficiency (absorption 
probability) is computed as the amount of radiation which is 
neither transmitted nor reflected. In terms of the elements of 
T, it is givenbya = l-(l + |Toi|2)|rn|-2 = 27(l + 7)-2. 

The curve shown in Fig. 2 reveals two regimes. If 7 <C 1, 
the decay channel |1) \g) is very slow compared to the 
time required for a photon to excite the |0) — > |1) transition, 
and only a small fraction of the photons is absorbed. If, on the 
other hand, the metastable state |1) decays too fast, 7^1, 
there is a Zeno suppression of the absorption. From previous 
formula, the maximal achievable detection efficiency is 50%, 
a limit reached by tuning the single absorber on resonance 
with the microwave field. We conjecture that this may be a 
fundamental limit for any setup involving a single point-like 
absorber and no time-dependent external control. 

A natural expectation would be that clustering many ab- 
sorbers inside the waveguide increases the detection effi- 
ciency. As shown in Fig. |2] this is not true. If we have a clus- 
ter of N identical absorbers close together, we can compute 
the detection efficiency using the same formula but with the 
scattering matrix Tdustor = . As far as the cluster size is 
smaller than a wavelength, the setup will be limited to a 50% 
maximum efficiency. There is a simple explanation for this. 
Since the cluster size is small, the photon sees the group of ab- 
sorbers as a single element with a larger decay rate, NT. This 
renormalization just shifts the location of the optimal working 
point, leaving the maximum efficiency untouched. 

The main result is that we can indeed increase the absorp- 
tion efficiency by separating the absorbing elements a fixed 
distance d longitudinally along the waveguide. The total scat- 
tering matrix for the array is given by 



N 

array — 6 

i=i 



d/ Xrp 



(6) 



where cr^ is a Pauli matrix and the scattering of each absorber 
may be different. In this case the microwave pulse does no 
longer see the the detection array as a big particle, and we 
obtain an collective enhacement of the absorption probabil- 
ity. Remarkably, an arbitrarily high detection efficiency can 
be reached by increasing the number of absorbers and tuning 
their separation d. Already with two and three qubits we can 
achieve 80% and 90% detection efficiency, see Figs. |2] and |3] 
Furthermore, as we have seen numerically, the more qubits we 
have, the less sensitive the whole setup becomes to the exper- 
imental parameters, see Fig. [2] This shows that the proposed 
setup is both robust and scalable. 

The previous analysis is of a general kind. It only requires 
a coupling between a waveguide and absorbers that may cap- 
ture a photon and irreversibly decay to one or more stable 
states. A practical implementation of our setup, which does 
not require strong coupling or cavities, consists on a copla- 
nar coaxial waveguide coupled to a number of current-biased 
Josephson junctions (CBJJ) 

HQ. We will now sketch the 
microscopic derivation of Eq. (O for this setup and relate the 
efficiency to the parameters of the circuit. 




FIG. 3: (a) Maximal achievable single-photon detection efficiency as 
a function of the number of qubits along the microwave guide, (b) 
Optimal working parameters vs. number of qubits. 



First of all, since the CBJJs are described by a washboard 
potential for the phase degree of freedom, we can identify |0) 
and |1) with the two lowest metastable levels in a local min- 
imum, see Fig.[T] The energy levels around such a minimum 
are well described by a harmonic oscillator with a frequency 
ijj that depends on the bias current. Furthermore, these levels 
have finite lifetimes before they decay into the continuum of 
current states, but since that the decay rate of state |0), Fq, can 
be made 1000 times smaller 11411 than that of Fi, we will 
approximate Fq ~ 0, Fi = F. 

The microwave guide is described by a Lagrangian |0] 



L = I dx 



(7) 



where / and c are the inductance and capacitances per unit 
length. The quantization of the charge field Q introduces 
Fock operators associated to the normal modes of the line. 
If we assume periodic boundary conditions, then Wp{x, t) ~ 
cxp[i{px — ujpt)]/ L~^^'^ , where L is the length of the waveg- 
uide and the dispersion relation is ujp = Vg\p\ = \p\/\/ci. 
When the relevant modes of the electromagnetic field are con- 
centrated in a small interval B around a central momentum 
Pq, we can introduce right and left moving fields tprix, t) = 
Epee°^p'"p(^'*)' ^ndijji{x,t) = Y,^^^a-pW-p[x,t), and 



pUjpa^ap, 



approximate the waveguide Hamiltonian as H 
which corresponds to the first line in Eq. Q. 

Finally, for the interaction between the absorbers and the 
guide we use a capacitive coupling in the dipole approxima- 
tion iSl . The corresponding Hamiltonian has the form 



Co 



C„ + Cj 



./^(^^+^,+H.c.)-(a + at), (8) 
V c a 



The first fraction depends on the capacitances of the junction 
and of the gate between the junction and the microwave, Cj 



4 



and Cg, respectively. The second term gives the strength of 
the electric potential inside the waveguide and it is propor- 
tional to the fields. Finally, the third term is just the charge 
operator for the CBJJ expressed using harmonic approxima- 
tion around a minimum of the washboard potential. In par- 
ticular, a ~ |0) (1| and = Ec/huj is the dimensionless 
parameter of this oscillator, expressed in terms of the junc- 
tion capacitance, Eq = (2e)^/Cj, and the plasma frequency. 
Note that when we combine all constants to form the interac- 
tion strength V there is not explicit dependence on the length 
of the microwave guide. Qualitatively, while in cavity experi- 
ment the qubit only sees a small fraction of the standing waves 
with which it interacts, in our setup each absorber gets to see 
the whole of the photon wavepacket after a long enough time. 

In terms of the microscopic model, it is possible to write the 
parameter that determines the detector efficiency as follows 



h Ti — i{u! - 



ci2 eZo 



(9) 



where we have introduced the dimensionless constant ci2 = 
Cg/ {Cg + Cj). It is evident from Eq. (|9]l that, in order to op- 
timize the efficiency, we have several experimental knobs to 
play with. In particular we have considered the following val- 
ues, close to current experiments ||20| a junction capacitance 
of Cj = 4.8pF, ci2 = 0.13 and w = 5 GHz. Putting the 
numbers together, and letting the waveguide impedance oscil- 
late between 10 and 100 SI, the optimal operation point for a 
single junction gives a necessary decay rate F ~ 10 — 100 
MHz. Increasing Cg by a factor 2 triples the optimal decay 
rate, F ~ 30 - 300 MHz. 

Our proposal has the following potential hmitations and im- 
perfections. First, the bandwidth of the detected photons has 
to be small compared to the time required to absorb a photon, 
roughly proportional to 1/F. Second, the efficiency might be 
limited by errors in the discrimination of the state \g) but these 
effects are currently negligible |25-| . Third, dark counts due to 
the decay of the state |0) can be corrected by calibrating Fo 
and postprocessing the measurement statistics. Fourth, fluctu- 
ations in the relative energies of states |0) and |1) , also called 
dephasing, are mathematically equivalent to an enlargement 
of the incoming signal bandwidth by a few megahertz and 
should be taken into account in the choice of parameters. Fi- 
nally and most important, unknown many-body effects cause 
the non-radiative decay process 1^0, which may manifest in 
the loss of photons while they are being absorbed. In current 
experiments 1251 this happens with a rate of a few megahertzs, 
so that it would only affect very long wavepackets. 

A natural extension of our design is the implementation of 
a photon counter. This only requires a number of absorbers 
large enough to capture all the photons that have to be detected 
and counted. Furthermore, while we have chosen three-level 
systems and CBJJs for illustrative purposes, the same goal 



may be achieve using other level schemes and quantum cir- 
cuits that can absorb photons and irreversible decay into long 
lived and easily detectable states. 
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